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ABSTRACT 

We report on our Spitzer observations of the anomalous X-ray pulsar 4U 0142+61, made following 
a large X-ray burst that occurred on 2007 February 7. To search for mid-infrared flux variations, four 
imaging observations were carried out at 4.5 and 8.0 /jm with the Infrared Array Camera from February 
14 to 21. No significant flux variations were detected, and the average fluxes were 32.1+2.0 /iJy at 
4.5 /im and 59.8+8.5 /iJy at 8.0 /im, consistent with those obtained in 2005. The non-detection of 
variability is interesting in light of reported rapid variability from this source in the near-infrared, but 
consistent with the fact that the source already went back to its quiescent state before our observations 
began, as indicated by contemporaneous X-ray observations. In order to understand the origin of the 
near-infrared variability, frequent, simultaneous multi-wavelength observations are needed. 

Subject headings: infrared: stars — pulsars: individual (4U 0142+61) — X-rays: stars — stars: 
neutron 
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1. INTRODUCTION 

The Anomalous X-ray Pulsars (AXPs) are a small 
group of isolated, young neutron stars having X-ray 
luminosities greater than their rotational energy loss 
rates. They are believed to be magnetars, neutron stars 
possessing ^10 14 G magnetic fields, with X-ray emis- 
sion powere d by the decay of their su perstrong mag- 
netic fields (jThompson fc Duncan! [1996T ). The magnetar 
nature of AXPs has been strongly supported by their 
high spin-down rates and similar radiative properties 
to soft gamma repeaters, i ncludin g sho rt X-ray bursts 
(jWoods fc Thompsonll200l lKaspill2007h . The bursting 
activity likely reflects structural changes in the surface 
magnetic field of a magnetar. 

As the brightest and nearest among the currently 
known AXPs, 4U 0142+61 has been well studied over 
wavelength ranges from the X-ray to mid- infrared (MIR) . 
Its X-ray emission in the 0.5-10 keV range can be well 
described by a blackbody (kT = 0.5 keV) plus po wer law 
(photon index T = 3.4; e.g., iPatel et al.ll2003f) . which 
is believed to arise fro m the surface and magnetosphere 
of t he star (e.g., iThompson. Lvutikov. fe Kulkarnil 
2002). Its opti cal emission is pulse d at the 8.7- 
s sp in period (jKern fc Martini l2002t iDhillon et all 
l2005h and appears to have a power-law-like spec- 
trum (Hullcma n. van Kerkwijk. fc K ulkarni 2000), likely 
also originating from the magnetosphere (although no 
apparent connection can be established between the 
optical and X-ray spectra). A surprise came from 
the detection of the source in the MIR (4.5 and 
8.0 fan), which reveale d a rising spectrum from the 
near- IR (NIR) to MIR (jWang. Chakrabartv. fc Kaplan! 
2006). The detectio n can be in t erpre ted as a sur- 
rounding debris disk (|Wang et al.l I2006D . The puta- 
tive disk, presumably formed from fallback material af- 
ter the supernova explosion, is irradiated by the X- 
rays from the central pulsar and emits mainly in the 
MIR. The existence and appearance of such a disk 
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has been predicted (e.g., iLin. Wooslev. fc BodenheimeH 
119911 : iPerna. Hernquist fc Naravanll2000l ). 

Recently, optical/NIR fl ux variability from 
4U 142+61 has been reported (|Durant fc van Kerkwiikl 
2006) . In a total of only 9 observations, significant flux 
variations were detected, suggesting such variability 
is very common. This variability is intriguing, since 
no such flux changes have been detected in X-rays, 
even though X-ray observations of the source have 
been made much more frequently than the optical/NIR 
observations. An obvious question is whether there 
is similar variability in the MIR. Rapid variability in 
the MIR would be very hard to understand in the disk 
model in the absence of X-ray variability. 

On 2007 February 7, a large fast-rise X-ray burst 
was detected during the Rossi X-Ray Timing Ex- 
plorer (RXTE) mon itoring observations of 4U 0142+61 
i|Gavriil et"all I2007T ). The peak flux of the burst was 
more than two orders of magnitude larger than that in 
the source's quiescent state. This event provided an op- 
portunity to test the fallback disk model, since a MIR 
flux increase is expected from an X-ray irradiated disk 
when the input X-ray flux is increased (see Figure [lj . 
For the purposes of searching for variability and testing 
the fallback disk model, we observed 4U 0142+61 with 
the Spitzer Space Telescope. In this paper, we report on 
the results of the observations. 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. Spitzer/IRAC 4-5/8.0 /im Imaging 

We observed 4U 0142+61 four times in 2007 February 
14-21 with Spitzer, following the February 7 burst. To 
catch flux variations from the source and possibly con- 
strain their time scale, the observations were scheduled to 
be on days 2, 4, and 7 after the first observation. The ex- 
act observation dates are given in Table [TJ The imaging 
instrument used was the Infrared Array Camera (IRAC; 
iFazio et al.1 l2004h . It operates in four channels at 3.6, 
4.5, 5.8, and 8.0 /im, while two adjacent fields are simul- 
taneously imaged in pairs (3.6 and 5.8 /zm; 4.5 and 8.0 
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/«ru). We observed our target in the 4.5 (bandwidth 1.0 
/urn) and 8.0 fim (bandwidth 2.9 /zm) channels. The de- 
tectors at the short and long wavelength bands are InSb 
and Si: As devices, respectively, with 256x256 pixels and 
a plate scale of T/2. The field of view (FOV) is 5.'2x5.'2. 
The frame time was 100 s, with 96.8 and 93.6 s effective 
exposure time per frame for the 4.5 and 8.0 /Ltm data, re- 
spectively. The total exposure times in each observation 
were 19.4 min at 4.5 /im and 18.7 min at 8.0 /jm. 

We also obtained the previous Spitzer IRAC data of 
4U 0142+61 from the Spitzer archive. The observation 
was made on 2005 January 17 at the same wavelength 
bands as ours. The effective exposure times were 77.4 
and 74.9 min in the 4.5 and 8.0 /zm, respectively. The 
results from this ob servation were previously reported by 
I Wang et all (p006h . 

2.2. Data Analysis 

The raw image data were processed through the IRAC 
data pipelines (version S15.3.0) at the Spitzer Science 
Center (SSC). In the Basic Calibrated Data (BCD) 
pipeline, standard imaging data reductions, such as re- 
moval of the electronic bias, dark sky subtraction, flat- 
fielding, and linearization, are performed and individual 
flux-calibrated BCD frames are produced. In the post- 
BCD (PBCD) pipeline, radiation hits in BCD images are 
detected and excluded, and BCD frames are then com- 
bined into final PBCD mosaics. The detailed data re- 
ductions in the pipeline s can be found in t he IRAC Data 
Handbook (version 3.0; lReach et alj |2006). For our 2007 
data, the target's nearby region in the 4.5 fim images suf- 
fered the "column pull-down" effect, an artifact in which 
a column of the detector array has a reduced intensity 
when a place in the column reaches a level of approxi- 
mately 35000 digital numbers (DNs) due to a bright star 
or cosmic rays. We corrected the BCD images for the 
artifact by using the tool cosmetic.pl. The artifact- 
corrected BCD images were then combined into PBCD 
mosaics by using the tool mosaic . pi. Both tools are pro- 
vided by SSC in the M0PEX package. In our PBCD data 
processing, the latest June 2006 permanently damaged 
pixels mask (pmask) image, provided by SSC, was used. 

The target field in the 4. 5 /xm PBCD i mage s is 
crowded. We used D0PH0T (jSchechter et all 11993( 1 . a 
point-spread function (PSF) fitting program, for pho- 
tometry. For the 8.0 /im images, in which the sky back- 
ground dominates in brightness and varies substantially 
over the source field region, we found that aperture pho- 
tometry provided more consistent measurements (this 
was also suggested by the Spitzer Helpdesk). The pro- 
gram phot in the IRAF package apphot was used for 
photometry. The aperture radius was 3 pixels and back- 
ground annulus radii were 3-7 pixels. Th e D0PH0T re- 
sults were also corrected to this aperture set (jReach et al.l 
2006). The correction factor was derived by comparing 
the brightnesses resulting from the PSF-fitting photom- 
etry to those from the aperture photometry of a few well 
exposed, isolated stars in the source field. 

In addition, to eliminate systematic variations among 
the PBCD images that might be caused by changes such 
as in the instrument, observing conditions, or data re- 
duction, we applied a differential photometry technique: 
we calibrated the brightness of our target with those of 
an ensemble of bright, non-variable stars in the images. 



The base images were the 2005 data, which have longer 
exposure times and thus higher signal-to-noise ratios. 




10 



20 



30 



40 



50 



/ in -10 -1 -2\ 

{ 1 ergs s cm ) 



Fig. 1. — The 4.5 fim flux (solid curve) as a function of the X-ray 
flux (0.5-10 keV) d erived from the dust disk model for 4U 0142+61 
llWang et all 120031 . Similar NIR K band (2.2 fim) flux changes 
(dashed curve) are also expected if the disk structure is assumed 
to be fixed. 

3. RESULTS 

The flux density measurements are given in Table [TJ 
As can be seen, the flux densities are consistent with be- 
ing a constant within the uncertainties. No significant 
variations were detected from our observations. Com- 
pared to the 2005 flux densities, the upper limits (90% 
confidence) on the flux density changes are 13% at 4.5 
/zm and 38% at 8.0 /im. We note that the flux densi- 
ties in Tabled] for the 2005 observation are slightly lower 
than those previously reported, 36.3 ±3.6 and 51.9±5. 2 
iijy at 4.5 and 8.0 (im, respectively (jWang et al.ll2006h . 
For the 4.5 /xm results, the difference is mainly caused by 
blending with nearby stars. A few photons from nearby 
stars appear within 2-3 pixels radii from the center of 
our target. Whe n aperture photometry was used (in 
I Wang et al.ll2006l ). the extra photons were included as 
the target counts. This can be verified by comparing the 
measurements from aperture photometry to those from 
PSF-fitting. For example, when an aperture radius of 
3 pixels is used, the ratio between the fluxes from the 
first and latter photometry is 1.42±0.09 for our target, 
while the average ratio is 1.214±0.008 for a few non- 
blended test stars in the field. When the aperture radius 
is changed to 2 pixels, the flux ratios for our target and 
the test stars are equal (1.07±0.07 for the target versus 
1.070±0.008 for the test stars). For the 8.0 /xm results, 
the difference, which is well within the uncertainties, is 
probably due to the differ e nt IR AC data pipelines [ver- 
sion S11.0 in IWang et all ((20061 ) versus S15.3.0 in this 
work] . 

We average the four 2007 4.5 fim fluxes, weighted by 
the flux uncertainties, and find a mean of 32.1±1.2 /iJy. 
This value is exactly the same as the 2005 result (note 
that the total exposure time in each wavelength band 
of the 2007 observations is approximately the same as 
that of the 2005 observation). Adding a 5% absolute 
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calibration uncert ainty in quadrature (|Reach et al.ll2005t 
iReach et al.ll2006t ). we obtain 32.1±2.0 /xJy. The 8.0 ^m 
flux densities, as shown in Table [TJ have similar uncer- 
tainties, even though the exposure time of the 2005 ob- 
servation is four times that of one 2007 observation. This 
is because the flux de nsity uncertainties are dominated 
by the sky noise (e.g. jReach et al.ll2005f ) and there were 
intrinsic intensity variations within the sky annulus in 
the aperture photometry. Therefore, the longer 2005 ex- 
posure time did not effectively reduce the sky noise. We 
average the four flux densities and find an average of 
59.8±8.5 fiJy, where we use the uncertainty 8.0 /iJy (Ta- 
ble QJ and add the 5% absolute calibration uncertainty 
to it in quadrature. 

4. DISCUSSION 

There are a total of 11 reported -ftT-band obser- 
vations of 4U 0142+61 in the past 9 years, and 
three of them were with simu l taneous JH observat ions 
(jDurant fe; van Kerk wiik 20061 [Gonzalez et 8,1.112007^ 1 . In 
the observations, the AXP exhibited large and rapid flux 
variations in the K band, and the JH fluxes correlated 
with the K but varied less. Assuming that the -ftT-band 
variations are sufficiently sampled, the probability of see- 
ing a >15% flux change would be 60%. In one extreme 
case, a 40% K flux decrease occured within one day. Our 
Spitzer observations were planned so that such strong 
variability would have been detected. Correlated flux 
changes between the MIR and K bands are expected in 
the proposed dust disk model for 4U 0142+61 (FigureHJ 
iWang et al.ll200l . if X -ray flux changes are considered 
as the primary cause of the variability. 

We did not detect the K-\ike flux variations in our 
MIR observations. The Spitzer /1RAC fluxes were con- 
sistent with being a constant within the uncertainties. 
Given that the upper limit on the MIR 4.5 (im flux vari- 
ations was < 13%, the MIR emission variability was sim- 
ilar to the NIR JH bands and less variable than in the 
K band. We note that in a i-T-band observation made 
one day prior to our first Spitzer observation, t he source 
had K — 19.9 ±0.1 mag (|Gonzalez et a,l.lf2007h . approx- 
imately the average K magnitude of 4U 0142+61. Be- 
cause the time scale of the K flux variations is not deter- 
mined, this could suggest that the source was in a stable 
state in the NIR during our observations, which could be 
a reason for the absence of MIR variability. To further 
constrain the relation between the MIR and NIR emis- 
sion, frequent, simultaneous observations are needed. 



If the MIR emission is truly less variable, the restricted 
-ft"-band variability would demand explanation. Interest- 
ingly, similar NIR variability is s een in a few dust disk 
systems around young stars (e.g.. lEiroa et al.ll2002l ). In 
those systems, the variability has been suggested to be 
caused by structural changes of the inner disk. If this is 
the case for 4U 0142+61, the properties in the source's IR 
variability might be explainable. Since emission in the K 
band would mainly arise from the inner edge of the disk 
and thus be more sensitive to inner disk changes, it would 
vary more strongly than emission in the JH and Spitzer 
MIR bands. It is not clear what would be the cause of 
the structural changes in the disk in 4U 0142+61. One 
possibility is variable X-ray heating, such as due to burst- 
related X-ray brightening. However, no significant short- 
term flux variations hav e been found in the X-ray ob - 
servations of the source (|Durant Sz van Kerkwiikl f2006). 
and thus far only a few bursts have been fo und in bi- 
week l y RXTE X-ray mon itoring observations ()Dib et al.l 
l2006t [Gavriil et ai1)2007h . indicating infrequency of the 
burst events. Thus in the absence of X-ray variability, 
the if-band variability challen ges the dust disk model 
(jDurant k, van Kerkwijkll2006h . 

In addition, we did not detect any significant MIR flux 
changes from 4U 0142+61 in our Spitzer IRAC obser- 
vations following the large X-ray burst. The 4.5/8.0 
/*m fluxes were nearly the same as those obtained in 
2005. As demonstrated by X-ray observations of the 
source that were made several days prior to our first 
Spitzer observation (Gonzalez et al. 2007), the AXP 
unfortunately had already returned to quiescence dur- 
ing our observations. Also as shown in the long-term 
X-ray monitoring observations, the quiescent X-ray flux 
of 4U 0142+61 h as been nearly constant and stable over 
the past decade ([Dib, Kaspi. fc Gavriill [2007; Gonzalez 
et al. 2007). Therefore, our non-detection is consistent 
with the source's state at X-ray energies at the time. 
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TABLE 1 

Spitzer IRAC 4.5/8.0 fiM flux density 

MEASUREMENTS OF AXP 4U 0142+61 



Observation Start Time 


F4.5 


^8.0 


(UTC) 


(**Jy) 


(**Jy) 


2005-01-17 22:31 


32.1±1.2 


48.8±7.6 


2007-02-14 05:50 


31.2±2.2 


58.2±8.8 


2007-02-15 17:59 


33.0±2.5 


53.5±8.8 


2007-02-18 09:02 


32.8±2.2 


65.6±8.0 


2007-02-21 15:58 


31.5±2.3 


61.2±9.8 



Note. — 5% absolute calibration uncertainty 
is not included. 



